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ABSTRACT 

We compare the angular expansion velocities, determined with VLBI, with the linear expansion 
velocities measured from optical spectra for supernova 1993 J in the galaxy M81, over the period from 
7 d to ^9 yr after shock breakout. The high degree of isotropy of the radio shell's expansion, within 
5.5%, with the projection of the radio shell being circular within even 1.4%, and the consistency of 
the radio shell thickness with predictions from hydrodynamic simulations and analytical computations 
argue strongly in favor of the radio-emitting shell being bound by the forward and reverse shocks. 
The absorption and emission of hydrogen at the highest velocity most likely also arise between the 
forward and reverse shocks, specifically in the shocked ejecta behind the contact surface which is 
close to the reverse shock, and possibly in Ray leigh- Taylor fingers that extend beyond it, farther 
out into the radio shell but not beyond the forward shock. The radio shell and the Ha absorbing 
and emitting gas are on average similarly decelerated, but the latter slightly less so than the former 
several years after shock breakout. This may indicate developing Ray leigh- Taylor fingers, extending 
progressively further into the shocked circumstellar medium. We estimate the distance to SN 1993J 
using the Expanding Shock Front Method (ESM). We find the best distance estimate is obtained by 
fitting the angular velocity of a point halfway between the contact surface and outer shock front to 
the maximum observed hydrogen gas velocity. We obtain a direct, geometric, distance estimate for 
M81 of D = 3.96 ± 0.05 ± 0.29 Mpc with statistical and systematic error contributions, respectively, 
corresponding to a total standard error of ±0.29 Mpc. The upper limit of 4.25 Mpc corresponds to 
the hydrogen gas with the highest observed velocity reaching no farther out than the contact surface a 
few days after shock breakout. The lower limit of 3.67 Mpc corresponds to this hydrogen gas reaching 
as far out as the forward shock for the whole period, which would mean that Rayleigh- Taylor fingers 
have grown to the forward shock already a few days after shock breakout. Our distance estimate is 
9 ± 13% larger than that of 3.63 ± 0.34 Mpc from the HST Key Project, which is near our lower limit 
but within the errors. 

Subject headings: supernovac: individual (SN 1993J) — radio continuum: supernovae 



1. INTRODUCTION 

Supernova 1993J was discovered in a spiral arm of M81 
south south-west of the galaxy's center by Garcia (Ripero 
& Garcia, 1993) on 28 March 1993, shortly after shock 
breakout at ^ UT (Wheeler et al. 1993) on the same 
day {t = d). It subsequently became the optically 
brightest supernova in the northern hemisphere since SN 
1954A and one of the brightest radio supernovae ever 
detected. It is also one of the closest extragalactic su- 
pernovae ever observed and is second only to SN 1987A 
as a subject of intense observational and theoretical su- 
pernova studies. This combination allows for sensitive 
spectroscopic and spectropolarimctric observations and 
for exceptionally detailed VLBI imaging, the latter giv- 
ing us the highest relative resolution ever obtained for 
any radio supernova (Figure [T]) . 



Combining the radial velocities of the ejecta gas ob- 
tained from the optical lines with the transverse angular 
velocities of the radio shell obtained from VLBI mea- 
surements yields a direct estimate of the distance to the 
supernova and its host galaxy. This method, called the 
Expanding Shock Front Method (ESM), was used to de- 
termine the distance to SN 1979C in MlOO in the Virgo 
cluster (Bartel 1985; Bartel et al. 1985; Bartel & Bieten- 
holz 2003, 2005). However the relatively large distance 
to Virgo has made it hard to resolve the supernova in 
sufficient detail, and an accurate distance estimate for 
Virgo is still pending. 

In contrast, M81, the host galaxy of SN 1993J, is much 
closer. Its distance was recently determined via HST 
observations of Cepheids to be 3.63 ± 0.34 Mpc (Freed- 
man et al. 1994, see also Freedman et al. 2001, HST 
Key Project) and 3.93±0.26 Mpc (Huterer, Sasselov, & 
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Schechter, 1995). Here we report on a detailed compar- 
ison between the angular velocities obtained from the 
VLBI measurements of the radio shell radius and the ra- 
dial velocities obtained from the Doppler shifts of Ha, 
H/3, He I, 0[ni], and Na I optical lines and determine 
the distance to SN 1993 J and its host galaxy M81 with 
ESM. 

This paper is the fourth in a series, presenting the re- 
sults from our VLBI campaign on this supernova (see 
Bartel et al. 1994 for early results and Bartel et al. 2000 
for an introduction of this series of results, see Marcaide 
et al. 1997 for parallel observations). In the first paper 
(Bietenholz, Bartel, & Rupen 2001; Paper I), we located 
the explosion center with respect to the nuclear radio 
source of M81, thus defining a stable reference point for 
our images. We also determined, using model-fitting, the 
motion of the geometric center of SN 1993J, and com- 
ment on the high degree of circular symmetry shown by 
SN 1993J. In the second paper (Bartel et al. 2002; Pa- 
per II), we determined the expansion speed of SN 1993 J 
and measured its deceleration as a function of time. In 
the third paper, we presented a complete series of VLBI 
images of SN 1993J at 8.4 and 5.0 GHz, along with our 
latest image at 1.7 GHz (Bietenholz, Bartel, & Rupen, 
2003; Paper HI). In this fourth paper we use the main re- 
sults of each of the previous papers, compare them with 
optical spectroscopic observations by others, include spe- 
cific predictions from hydrodynamic simulations of this 
supernova for comparison with our observations, and de- 
termine directly the distance to M81. 

In § [21 we summarize our VLBI observations and data 
reduction. In §|3]we describe the generic supernova shell 
model with its radio and optical emission regions and dis- 
cuss the origin of the radio and optical emission and the 
relation of the forward and reverse shock to the measured 
boundaries of the radio emitting shell. In §[4] we elabo- 
rate on the expansion of the supernova as measured in 
the radio and optical wavelength range. Then we deter- 
mine the distance to SN 1993 J and its host galaxy in §[5l 
discuss our results in § [S] and finally give our conclusions 
in m 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were described in Papers I and II. To 
summarize, we observed SN 1993J at 34 epochs between 
1993 and 2001 from i = 30 to - 3000 d. At the earliest 
epoch, at i = 30 d, we observed only at 22.2 GHz. At 
later epochs up to the 33'''^ epoch at i = 2996 d, we 
observed mostly at 8.4 GHz and often also at 5.0 GHz. 
In addition, we also observed at some early epochs at 14.8 
GHz and throughout the full eight years sporadically also 
at 2.3 and 1.7 GHz. We used a global array of between 
9 and 18 telescopes with a total time of 9 to 18 hours 
for each session. We have continued observing SN 1993J 
at approximately one-year intervals and will report on 
these observations in future papers. 

The data were recorded with either the MK III or the 
VLBA/MKIV VLBI systems, and correlated with the 
NRAO VLBA processor in Socorro, New Mexico, USA. 
We refer the interested reader to Paper II where details 
of the observing sessions are tabulated. The analysis was 
carried out using NRAO's Astronomical Image Process- 
ing System (AIPS). 

Most of our observations were phase-referenced to 



M81*, the radio source in the center of the galaxy M81. 
Being very compact (Bietenholz et al. 1996), M81* is an 
excellent calibrator for phase-referenced mapping. Our 
images are among those with the lowest background noise 
level currently obtained with VLBI. Bietenholz, Bartel, 
& Rupen (2000) located a fixed point in the variable 
brightness distribution of M81* that also has an inverted 
radio spectrum (Bietenholz, Bartel, & Rupen 2004a). 
We identified this point with the purported supermas- 
sive black hole in the center of the galaxy. Using this 
point as a reference point allowed us to locate the explo- 
sion center of SN 1993J on each of our supernova images 
with high precision. 

In Figure [1] we show the galaxy M81 at optical wave- 
lengths. Overlayed is an image of the galaxy at radio 
wavelengths, with SN 1993J clearly visible in a south- 
ern spiral arm. On the right side of the figure we show 
the VLBI image of the evolved radio shell which has the 
highest relative resolution obtained so far (Paper III). 

In addition to imaging, we used model-fitting. We fit, 
by weighted least-squares, the two-dimensional projec- 
tion of a three-dimensional spherical shell of uniform vol- 
ume emissivity to the calibrated u-v data consistently for 
all epochs. The ratio of the outer to inner angular radius 
was fixed at 9o/9i— 1.25. From this fit, we estimated the 
shell's center coordinates, x, y, (see Paper I), and 9o (see 
Paper II). The shell's center coordinates were found to be 
equal, within an rms of 64 /las, to those of the explosion 
center of SN 1993J. In the image in Figure [T] this center 
is at the origin. In addition, for relatively late epochs, we 
also determined the shell thickness by freeing the ratio 
of the radii and estimating 6o and 9i independently. 

3. THE GENERIC SUPERNOVA SHELL MODEL AND ITS 
RELEVANCE FOR ESM 

To clarify the astrophysical interpretation of the ra- 
dio shell and its relation to the optical absorption and 
emission by the ejecta, we copied the left part of the 
VLBI image from Figure [T] and juxtaposed it to a sketch 
of a generic supernova shell model in Figure [2l The 
model describes a supernova expanding into the circum- 
stellar medium (GSM) left over from the progenitor star. 
The supernova ejecta are spherically, freely expanding 
away from the explosion center for a large range of radii, 
r, from that center. Their temperature is 10^ K. 
The ejecta in the outermost regions hit the GSM, and 
a contact surface forms between the two media at a ra- 
dius, Tcs. From this surface, a forward shock, located 
at a radius, rfs, travels into the GSM and heats it up 
to ~ 3 X 10* K. Outside the forward shock, the GSM 
is unshocked and has a relatively high temperature of 
10^ to 10^ K (Lundqvist & Fransson 1988; Fransson 
& Bjornsson 1998; Mioduszewski, Dwarkadas, & Ball 
2001). Also, a reverse shock travels from the contact 
surface back into the ejecta to a radius, r^-g, from the 
explosion center and heats up the shocked ejecta to a 
temperature of ^ 10*" K. The contact surface itself is 
Rayleigh- Taylor unstable, and fingers of ejecta are ex- 
pected to develop and extend into the shocked GSM. 
These fingers are also sketched in Figured] 

This general scenario above was described mathemat- 
ically by the mini-shell model of self-similar expansion 
(Ghevalier 1982a, b; see also Fransson, Lundqvist, & 
Ghevalier 1996). This model assumes power-law density 
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profiles for both the ejecta and the CSM, with pcj (x r " 
{n >5), and pcSM « r"^. With this assumption, the ex- 
pansion is self-similar, and the shock radius is given by 
Tfs ex t™, where m is the deceleration parameter, which 
is constant and is given by m = (n — 3)/(ri — s). A 
hydrodynamic model for SN 1993J by Mioduszewski et 
al. (2001), based on an ejecta model by Shigeyama et 
al. (1994) (see also Suzuki & Nomoto 1995), relaxed the 
assumption of a power-law density profile in the ejecta. 
This model described the expansion of SN 1993J in more 
detail, and accounted for the deviations from self-similar 
expansion which were found for this supernova (Paper 
II). 

The Expanding Shock Front Method is based on the 
transverse velocity of the Rayleigh- Taylor unstable con- 
tact surface being equal to the largest radial velocity of 
the ejecta. Dividing the latter by the former gives a ge- 
ometric determination of the distance to the supernova 
and its host galaxy. In the following we will describe 
how the contact surface velocity distorted by Rayleigh- 
Taylor instabilities is related to the expansion velocity of 
the outer surface of the radio shell and how it compares 
with the maximum velocities measured in optical lines of 
the ejecta. 

3.1. The origin of the radio emission 
3.1.1. The forward and reverse shocks 

It is believed that the forward shock accelerates elec- 
trons to ultrarelativistic velocities. These relativistic 
electrons then interact with the magnetic field amplified 
near the Rayleigh- Taylor unstable contact surface (e.g.. 
Chevalier, Blondin, & Emmering 1992; Jun & Norman 
1996a, b). This interaction leads to radio synchrotron 
radiation. Radio emission is therefore expected to orig- 
inate in the shell region of the shocked CSM and likely 
also the shocked ejecta, i.e., in a region bound by r^g 
and Tfg. No radio emission is expected from the freely 
expanding, unshocked ejecta. Radio emission could also 
originate from a central compact source associated with 
the stellar remnant of the explosion, namely a neutron 
star or a black hole. Although a compact source was 
indeed found inside the shell of SN 1986J (Bietenholz, 
Bartel, & Rupen 2004b) , no emission from such a source 
has yet been detected for SN 1993J or any other modern 
supernova (e.g., Bartel & Bietenholz 2005). 

How can our radio images of SN 1993J be interpreted 
in detail in the context of this astrophysical model? The 
radio image of SN 1993 J in Figure [T] shows an exemplary 
shell of emission. In general, the outer angular radii, 9^ 
of our geometric shell model could be determined with 
an accuracy of about 1 to 2% and, for our latest images, 
the mean ratio of the outer to inner shell radii, Oo/Oi-, 
to within 1% (Paper II). In Figure [2] we plot the fit ge- 
ometric shell model as overlayed concentric circles with 
radii of 9o and 9i. The center of the model circles is at 
the explosion center within 64 /ias which corresponds to 
only about 1% of 0o (Paper I). Our finding that the radio 
emission can be so well fit by a spherical shell centered 
at the explosion center strongly supports the expectation 
that the fit angular radii, 9^ and 9i, indeed correspond 
to the radii of the forward and reverse shocks, rfs and 
r^s, respectively. Other supernovae may not show such 
exemplary shells of emission and would therefore not be 



as good examples to show the close relationship between 
the radio shell boundaries and the forward and reverse 
shocks although that relationship most likely also exists 
given the example of SN 1993 J. 

3.1.2. The thickness of the radio shell and its relation to the 
forward and reverse shocks 

Both the self-similar mini-shell model and the hydro- 
dynamic model make predictions as to the distance be- 
tween the forward and the reverse shocks. During the 
first months, when m = 0.92 (Paper II), the ratio be- 
tween the forward and reverse shock radii is predicted to 
be 1.19 for s = 2 for the self-similar mini-shell model. 
For a larger deceleration the shell thickness is predicted 
to be larger. For instance, with the mean of m for the 
whole observing time up to t ^ 3000 d of m = 0.83 (Pa- 
per II), a ratio of 1.30^ is predicted by the self-similar 
mini-shell model for s — 2 (Chevalier & Fransson 1994). 

Similar and most likely more reliable values for the 
shell thickness are predicted by the hydrodynamic model. 
In fact, Dwarkadas & Chevalier (1998) found that the 
shell thickness is relatively independent of the form of 
the ejecta profile for a long time after shock breakout, 
further supporting its robustness. In Figure [3] we plot 
the ratio of the forward and reverse shock radii from a 
few days to several years after the explosion as derived 
from the hydrodynamic simulations. The ratio increases 
from 1.16 between i = 10 d and 20 d to an average of 
1.29 between t = 1000 d and 3000 d. In addition we plot 
our measurements from VLBI observations, 9o/9i, deter- 
mined at epochs from t — 996 d onward when the radio 
shell was large enough to allow useful determinations. 
The mean of do/di for epochs up to i = 2996 d after the 
explosion was 1.29 ± 0.01, with a slight broadening of 
the shell (significant only at the la level, see Paper II), 
in good agreement with the models. This good agree- 
ment is additional evidence that the radio shell is indeed 
bound by the forward and reverse shocks. 

3.1.3. The brightness profile of the shell and its relation to 
the forward and reverse shock 

Below the image of the half radio shell in Figure [2] we 
plot the profile of the fit geometric shell model in several 
aspects. The data points give the observed profile. The 
solid curve gives the fit to the data of a (projected) shell 
model with uniform emissivity between 9i and 9o and 
partial absorption inside 9i (Paper III). The rectangular 
lines give the profile in the (unprojected) cross section 
of the shell. Clearly, the emission and the region where 
the emission originates is very well fit by our geometric 
model. In Paper III we reported on the steepness of the 
outer edge of the profile and found that it is indeed con- 
sistent with being an ideally sharp edge as in our model. 
However softer boundaries are not excluded. Neverthe- 
less, the measured brightness profile can be taken as fur- 
ther evidence that the radio shell is indeed bound by the 
forward and reverse shocks, and that the parameters 9i 
and 9o therefore most likely correspond closely to r^-g and 
ris- 

^ These are average values for SN 1993J. But since m is observed 
and s inferred to change with time, and the evolution is in general 
not self-similar (Paper II), the ratio may differ somewhat from the 
predicted value. 
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On the right side of Figure [2] we sketch the minisheU 
model and complete the concentric circles by indeed as- 
suming here and hereafter that 9i and 9^ correspond ex- 
actly to and rts. 

3.2. The origin of the optical emission 

While optical continuum emission of the supernova is 
expected to mostly originate, at least at early times, from 
the photosphere in the central region of the expanding 
gas, optical emission and absorption in broad spectral 
lines is expected to arise in the gas of the hot ejecta, 
above the cooler and mostly neutral ejecta (Figure [51 In 
particular, emission was discussed as coming from a) the 
freely expanding ejecta and being caused by the heating 
and ionizing radiation from the interaction region, as well 
as from b) the shocked ejecta and being linked to radia- 
tive cooling (Chevalier & Fransson, 1994, and references 
therein). The ionizing luminosity and the density of the 
gas both change with time and give rise to evolving spec- 
tra. Blending of lines is common, and since the relative 
luminosities change with time, the degree of blending can 
change as well. The highest expansion velocity in absorp- 
tion and emission lines is expected from the shocked gas 
of the envelope of the progenitor star, in a region that ex- 
tends outwards to the Rayleigh- Taylor unstable contact 
surface at radius res. The highest observed optical expan- 
sion velocity could be smaller. No optical line emission 
is expected from the region of shocked CSM, i.e., beyond 
the contact surface between res and rfs, since the tem- 
perature there is expected to be 3 x 10^ K. 

3.3. The contact surface 

3.3.1. The location of the contact surface within the radio 

shell 

The location of the contact surface in relation to the 
forward and reverse shocks, i.e., to the outer and inner ra- 
dio shell surfaces, is estimated by combining VLBI mea- 
surements of the radio shell thickness with circumstel- 
lar interaction models. Both the self-similar mini-shell 
model and the hydrodynamic model make predictions as 
to the distance between the contact surface and the re- 
verse shock. During the first months, when m = 0.92, 
the ratio between the contact surface and reverse shock 
radii is predicted to be 1.0065 for s = 2 by the self- 
similar mini-shell model. For a larger deceleration, e.g., 
for a mean oi m — 0.83, s = 2 as taken before for a com- 
parison fij I3.1.2|) . the ratio is 1.03 (Chevalier & Fransson 
1994). The hydrodynamic model predicts a ratio sim- 
ilarly close to unity during the first month which then 
increases to a mean value of about 1.04 between t ^ 1000 
and ~ 3000 d. 

We parametrized 9i for the range of our radio and 
the optical observations by 9i= Ai-ft™''' 6*0, where t is 
the time since explosion in days. With A^^ ~ 0.914 
and mrf = —0.0217, this parametrization gives a time- 
varying shell thickness consistent with the hydrodynamic 
simulations with, e.g., 9o/9i = 1.16 at i = 15 d and 1.29 
at t = 3000 d (see Figure [3]). These values are also fairly 
consistent with the predictions from the self-similar mini- 
shell model. (We also parametrized the angular radius of 
the contact surface by 9c= Acft"^"' 9o with Ad = 0.904 
and TOcf = —0.0150 to give 9c being 0.65% and 4% larger 
than 0i at t = 15 d and 2500 d, respectively, as predicted 
by the above models. 



3.3.2. Distortion of the contact surface by Rayleigh- Taylor 

instabilities 

However, since the contact surface is decelerating, 
it is expected to be Rayleigh- Taylor unstable (Gull 
1973; Chevalier 1982b), and that with time fingers of 
shocked ejecta will extend into the shocked CSM. A 
linear analysis of self-similar solutions as well as two- 
dimensional numeric hydrodynamic computations have 
shown that Rayleigh- Taylor instabilities are limited by 
Kelvin-Helmholtz instabilities but can still grow in length 
to generally 30 to 50% of the width of the shocked shell 
(Chevalier et al. 1992; Dwarkadas 2000). Under partic- 
ular circumstances, when vortices are amplified by the 
magnetic field, it was shown that the fingers may even 
reach slightly beyond the outer shock, which is otherwise 
assumed to be spherically symmetrical (Jun & Norman 
1996a, b). Also, for fast shocks with efficient particle 
acceleration and resulting high compression ratios, the 
width of the shocked shell was found to shrink consider- 
ably so that fingers, without growing much larger, could 
also reach and even slightly exceed the average forward 
shock radius (Blondin & Ellison 2001). 

3.3.3. The fastest parts of the ejecta 

The ejecta with the highest velocity are therefore lo- 
cated in the evolving Rayleigh- Taylor fingers stretching 
out beyond res toward the forward shock at rfg. The fin- 
gers are expected to be less developed a few days after 
shock breakout and more developed a few years there- 
after. Thus, the highest possible velocities in optical lines 
are those from the front of the Rayleigh- Taylor fingers^ 
which, depending on evolution time, stretch out to be- 
tween Tcs and rfs with an assumed mean, rx = ^ (rfs -|- res), 
which corresponds to an angular extent of 9^. 

4. THE EXPANSION 

4.1. The angular expansion velocities of different 
surfaces of and within the radio shell 

We used the values of 9^ obtained from our VLBI ob- 
servations to compute the angular expansion velocity of 
the outer surface of the radio shell, i.e., of the forward 
shock, 0o, as a function of time. The values of 9o were 
already listed and plotted in Paper II, but for the conve- 
nience of the reader we repeat the plot of 9o as a function 
of time here in Figure [H It can be clearly seen that the 
outer surface of the radio shell is only slightly deceler- 
ated up to i 300 d and from then on more strongly 
decelerated. Further, as can be seen from the inset of 
Figure m the deceleration changes significantly even af- 
ter t ^ 2000 d (see also Paper II). 

Because the expansion has an approximate power-law 
form, a linear fit to the values of 9o over a time interval 
would produce a biased estimate of 9^ at the midpoint 
of the interval. To avoid this bias, we computed 9o{t) by 
first fitting a running solution of the form 9o = At™ to 
the data points near t, where we fit both the decelera- 
tion parameter, m, and the angular size scale, A. Our 
running solution is equivalent to a running mean, but 
instead of computing the weighted average in a shifting 

^ Note, that whether such line emission is observable depends 
on the physical condition in the fingers. 
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interval we computed the weighted hnear fit in a shift- 
ing interval in a log-log diagram (see Paper II for more 
details). As our value of 9o{t), we take the derivative of 
this fitted function, Amt"^~^, evaluated at the midpoint 
of the fitting interval on a logarithmic scale"^. 

We extrapolated the angular velocity values slightly 
beyond the mean times of the first and last running so- 
lutions to cover the time range from t ~ 5 to ~ 3000 d. 
We plot the values of 6o and connect them by lines in 
Figure [Sj We also estimate standard errors from the 
power-law fits. These errors are between 1 and 5% of 6o- 
We plot them also in Figure [51 We further compute the 
angular velocity of the inner surface of the radio shell, i.e. 
of the reverse shock, Oi, the velocity of the contact sur- 
face, 9c, and the mean between the velocity of the contact 
surface and that of the forward shock, O^, and also plot 
them in Figure O The ratios of the parametrized veloc- 
ities, do/di and 9o/6c, are comparable to those from the 
hydrodynamic simulations after the latter are smoothed 
over appropriate time intervals.^ 

4.2. The radial expansion velocities of the optical 
line- emitting gas 

4.2.1. The data 

We list the maximum expansion velocities found in op- 
tical lines of hydrogen, oxygen, and sodium in Table [TJ 
We obtained the earliest velocities from measurements by 
Finn et al. (1995) who reported 17,800 declining to 16,600 
km s~^ from t = 7 to 10 d after shock breakout for the 
minimum of the Ha A6563 absorption trough. The max- 
imum velocities are larger. Since an absorption profile 
was not published, we needed to calculate the maximum 
velocities on the basis of absorption profiles at overlap- 
ping times pubhshed by others. For t = 10 to 19 d, maxi- 
mum velocities of 18, 000 km s^^ were obtained from the 
blue edge of the Ha absorption trough (Trammell et al. 
1993; Lewis et al. 1994). Since at i = 10 d the maxi- 
mum velocity is 1400 km s^^ higher than the velocity at 
the absorption minimum, we added 1400 km to Finn 
et al.'s values to obtain the maximum velocities also for 
t = 7 to 9 d. The highest velocity is 19,200 km s^^ at 
t = 7 d. 

At later times the absorption trough first shrank, most 
likely due to blending at the blue side with the rising 
[01] AA 6300, 6364 emission doublet^. Then the trough 
disappeared due to Ha becoming optically thin (see, e.g., 

^ In particular, we computed a "running solution," 9 = Ajvi™", 
for each N, from the values of do for 8 consecutive observing epochs 
from tjv to tjv+7. We started with the data from the first epoch 
at ti to the eighth epoch at tg for the first solution (A'' = 1), 
and continued till tjqj^-j was the last epoch. For each N, we fit 
for the values of yljv and mjv- The time at which these fitted 
values of Ajv and rrijv are taken to apply, tjy, is the geometric 
mean of the start and end times of the segment, i.e., logtjv = 
^(logtiv + logtjv+r)- At tjv we take the angular velocity to be 
given by do{t~N) = A^nifitM"'" ''^ . 

The ratio, rfg/rrs, predicted by the hydrodynamic simulations 
(Figure [Sjl shows an oscillation at late times which corresponds 
to an oscillation of the ratio rfs/rrs- However, our data are also 
consistent with our simple parametrized function which does not 
oscillate. Our ratios of the parametrized velocities at late times are 
similar to the mean of the high and low values of the oscillation of 
the velocity ratios from the hydrodynamic simulation. 

^ Patat, Chugai, &; Mazzali (1995) report a sharp decrease of 
the radial velocity of the absorption minimum in the Ha profile 



Houck & Fransson 1996), and the maximum velocity for 
Ha could only be derived from the Ha emission lines. 

Several authors have measured emission lines in 
SN 1993J, e.g., Lewis et al. 1994; Spyromilio 1994; Fihp- 
penko, Matheson, & Barth 1994; Clocchiatti et al. 1995; 
Patat et al. 1995; Houck and Fransson 1996; Matheson 
et al. 2000a, b; Fransson et al. 2005. Of these, Patat et 
al. (1995) and Matheson et al. (2000a, b) reported mea- 
surements over the longest time intervals. Patat et al. 
(1995) plot representative Ha profiles and list RVZI (red 
velocity at zero intensity) values for all epochs of their 
observations. We list these values in Table [1] Matheson 
et al. (2000a) plot the profiles for all epochs of their ob- 
servations and list (Matheson et al. 2000b) BVZI (blue 
velocity at zero intensity) as well as RVZI values. We 
considered the possibility of asymmetries of the line pro- 
files that could indicate absorption effects or possible bi- 
ases in the determination of the maximum velocities and 
therefore determined from the profiles consistently the 
BVZI and RVZI values^ as well as the BVHI (blue ve- 
locity at half intensity) and RVHI (red velocity at half 
intensity) values^. These values are also listed in Table[TJ 

4.2.2. The maximum Hot velocities 

For a supernova like SN 1993J, hydrogen gas is a domi- 
nant constituent of the outer regions of the freely expand- 
ing ejecta and the shocked ejecta. We therefore plot the 
BVZI, BVHI, RVZI, and RVHI values of the Ha profile 
from Table [T] as a function of time in Figure [6l 

It is apparent that the maximum Ha velocities decrease 
with time. In particular, the BVZI decreases at a rate 
similar to that of 9^ with the BVZI value a,t t = 2454 d 
being about half as large as that at t = 7 d. The BVZI 
values at early times, between t = 7 and 19 d, are deter- 
mined from the blue edge of the absorption trough. They 
are between 19,200 and 18,000 km s~^, the largest veloc- 
ities measured for SN 1993J at any time. After that time, 
the 0[I] A6300, 6364 line appears and steadily grows in 
intensity, rendering impossible a measurement of the Ha 
BVZI in either absorption or emission or of the BVHI in 
emission. Only at t > 523 d could the velocities on the 
blue side of the Ha line be measured again. Figure[n]fur- 
ther shows that the earliest RVZI values, from t 16 to 
~ 19 d, are ~20% smaller than the corresponding BVZI 
values. At later times, however, the RVZI values agree 
within the errors with the BVZI values. Furthermore, a 
relatively large discrepancy exists between the RVZI and 

from ~ 13,500 km s~^ to ~ 9500 km between t = 15 d to 
~ 50 d. From the spectra in this time range (Matheson et al. 
2000a) it appears that the sharp decrease is caused by the [OI] 
AA 6300, 6364 emission doublet filling in the blue side of the Ha 
trough even at these early times. 

® We calculate the velocity at zero intensity by extrapolating the 
side of the line through the noise to the zero intensity baseline of 
the spectrum. 

Matheson ct al. (2000b) list Ha BVZI and RVZI values but not 
BVHI and RVHI values. For consistency we did not use Matheson 
et al.'s BVZI and RVZI values but rather determined these values 
and the corresponding BVHI and RVHI values from the line pro- 
files given by Matheson et al. (2000a). Matheson et al.'s (2000b) 
values do not differ from our determined values significantly, ex- 
cept at t = 523 and 553 d where we think that their BVZI values 
are underestimates. These underestimates are likely caused by a 
misinterpretation of the baseline level, which we believe is taken 
too high due to the still present blend of the 0[I] A A 6300, 6364 
emission doublet (see §|A|. 
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RVHI values at early times, the latter being almost 60% 
smaller than the former. This discrepancy decreases at 
later times due to a steepening of the red side of the Ha 
line profile. At t > 500 d, the RVZI and RVHI values 
are different by only ~20% and match the corresponding 
values on the blue side, indicating a largely symmetric 
profile with steep red and blue sides. 

These characteristics of the Ha profile at times up to 
t ~ 500 d indicate that the red side may be biased by 
absorption, likely by dust mixed in with the ejecta. The 
occurrence of dust in supernovae and absorption by it 
was discussed by, e.g., Fransson et al. (2005), Gerardy et 
al. (2000), and Deneault, Clayton, & Heger (2003). The 
maximum velocity of the Ha line emitting gas is there- 
fore best revealed by the BVZI (absorption and emission) 
values at all times and the RVZI (emission) values for 
t > 500 d«. 

4.3. Radio versus optical deceleration 

How well does the deceleration of the radio shell match 
that derived from the maximum velocities from the opti- 
cal lines? In Table [2] we compare the deceleration param- 
eter of the radio shell, m, with that of the line absorbing 
or emitting gas, mopt, for the lines discussed above and 
for different time ranges. In particular we compute mo 
from 6*0 oc t™° , nic from 6*0 oc f"'', from O^oc f^^, and 
mopt from Vopt oc moptt'-™"'"'"^'' , where ■^opt is the maxi- 
mum velocity of the line absorbing or emitting gas. 

We first computed the deceleration of the radio shell 
and that of the Ha line absorbing and emitting gas for 
the total time range. Here we had to consider the un- 
even sampling of the radio and optical data. Since the 
deceleration is changing with time, a fit with a single 
deceleration parameter would depend strongly on the 
weighting scheme. In fact, on a logarithmic scale, the 
radio data for instance were only sparsely sampled at 
early times and more densely sampled at later times. If 
the non-uniform sampling is ignored, then as reported in 
Paper II, mo = 0.827 ± 0.004 for a fit weighted only with 
the data uncertainties as given in Table II of Paper II. 
However, such a fit represents a strong mismatch to the 
earliest data. For a better comparison with the optical 
data including the earliest data from t = 7 to 19 d we 
increased the weighting of the early radio data by forcing 
the fit through the earliest radio data point at i = 30 d. 
For the optical data such weighting was not necessary 
mostly due to the large gap at intermediate times. 

The resulting "average" deceleration parameter for the 
Ha gas is equal within the errors combined in quadra- 
ture to that of the outer surface of the radio shell and 
that of the surface between the latter and the contact 

* Patat et al. (1995) discussed a possible blending of the red side 
of the Ho emission profile at early times (4=255 d) with an uniden- 
tified feature reported to be present in Type Ib/c supernovae but 
limited the emission to not more than 30% of that of Ha. How- 
ever, this feature, visible in Patat et al.'s (1995) example, cannot 
be seen in the spectrum of SN 1993J. Further, Chevalier & Frans- 
son (1994) computed that another line, [N H] A6548-6583, could 
be blended with the Ha line. However, at t = 2 yr its luminosity 
would be only ^5% of that of the Ha line and at t = 5 yr only 
20 to 25%. These limits and the time frame of the occurrences 
of these possible blends together with the near consistency of the 
RVZI with the BVZI values for t > 500 d are strong indications 
that our values are not significantly affected by blending with such 
lines. 



surface (solution 1 in Table [2]). The equality within the 
errors also holds for the early time interval (solution 2 
in Table [5]). For the late time interval, the Ha line ab- 
sorbing and/or emitting gas is slightly, but significantly, 
less decelerated than the radio shell with a difference in 
the deceleration parameters of 0.092 (4(7) and 0.099 (4cr) 
between the optical line emitting gas on the one hand 
and the outer surface of the radio shell and the surface 
between the later and the contact surface, respectively, 
on the other hand (solution 3 in Table [5]) . 

Somewhat larger discrepancies in the sense of less de- 
celeration are found for the 0[III] and Na I line emit- 
ting gas. The largest discrepancy is found for the H/3 
line-emitting gas, where no significant deceleration was 
measured. This is due to the relatively large spike of the 
H/3 velocities at t ~2000 d (see Table [J). 

5. THE DETERMINATION OF THE GEOMETRIC 
DISTANCE TO SN 1993J AND M81 

5.1. The distance solution with statistical errors only. 

As we have argued, the most reliable measurements of 
the maximum velocity from optical lines spanning the 
longest time are those from the blue edge of the Ha ab- 
sorption trough at early times of t = 7 to 19 d and of 
the Ha emission profile later on. We therefore take these 
Ha BVZI values (Table [ij and fit to them 9^, which is 
the mean between 9c and 9o, taken to account for the 
Ray leigh- Taylor fingers. More precisely, we interpolated 
9o{t) to the times, tj, of the optical data, fopt(ij), and 
computed the corresponding O^itj) values to get distance 
estimates, 

^ ^ ^opt(^j) 

for each tj . We then took the weighted mean of the De- 
values, and its uncertainty, as the solution (solution 1 in 
Table [3 for the distance, D, to SN 1993J and therefore 
to M8F. 

£' = 3.96 ±0.05 Mpc. 

We plot the Ha emission and absorption BVZI values 
and the radio shell velocity curves including the velocity 
curve for the contact surface and the mean between it and 
the outer surface (see Figure [5|) for the newly determined 
distance of 3.96 Mpc in Figure [T] 

5.2. Sensitivity study 

To investigate how much our distance determination 
depends on our choice of the optical velocity measure- 
ments, we solved for D by using subsets of the Ha ve- 
locities, velocities from other lines, and combinations of 
velocity sets, all from Table [TJ We list the solutions also 
in Tabled 

^ If we ignore that the deceleration parameters for the radio 
shell and the hydrogen gas are changing with time, and take in- 
stead the "average" deceleration parameters, nix and rriopt, from 
solution 1 in Table (2] the distance is given as D = that 

is D = ^I^t^opt-^x with ,;„pt = Aoptmoptt'""?'-^ and 
6x= Axt"^^ . We get slightly lower values for D between 3.83 Mpc 
for i = 10 d and 3.91 Mpc for t = 2000 d. This bias is caused by 
the inferior fit to the radio data and the resulting bias of slightly 
higher predicted velocities at early and late times where the optical 
data were obtained. We do not consider these solutions further. 



Distance to M81 



7 



5.2.1. Distance with Ha velocities 

In particular, we used only the velocities from the blue 
edge of the Ha absorption profile, BVZI abs. (solution 
2 in Table [2]) , from the blue edge of the Ha emission 
profile, BVZI em. (solution 3), and from the red edge of 
the Ha emission profile, RVZI, from t — 523 to 2454 d 
(solution 4). These velocity values were already plotted 
in Figures El and [71 These distance estimates straddle the 
one in solution 1, differing from it by not more than 0.05 
Mpc or 1.3% (0.6cr, combined statistical uncertainty^"), 
having somewhat larger errors, and xt values closer to 
unity apart from xt = 0.20 in solution 2. 

5.2.2. Distance with H/S velocities 

We further used the velocities from the blue edge of 
the H/3 A4861 emission profile (BVZI) from t = 553 to 
2454 d (solution 5). We plot the H/3 BVZI values and for 
comparison, the Ha BVZI values with the radio velocity 
curves for the distance of 3.96 Mpc (solution 1) in Fig- 
ureH The H/3 BVZI values from t - 500 to - 1000 d are 
~10% smaller than the corresponding Ha BVZI values, 
however they get larger with time and a± t ^ 2000 d, 
even exceed the Ha BVZI values, just at the time when 
the outer radio shell velocity is spiking up. It is inter- 
esting to note that this behavior is reflected in a bump 
in the lightcurve at X-rays (Zimmermann & Aschenbach 
2003). The distance with these values only (solution 5) 
is within 0.06 Mpc or 1.5% (0.3(7, combined statistical 
uncertainty) equal to the distance from solution 1. 

5.2.3. Distance with all hydrogen velocities combined 

Since the hydrogen gas constitutes a large fraction of 
the outer parts of the ejecta, we combined all the BVZI 
and RVZI values from the Ha and H/3 profiles used for so- 
lutions 1, 4, and 5 and again solved for the distance. The 
resulting distance (solution 6) is only 0.01 Mpc or 0.3% 
(0.2ct, statistical uncertainty from solution 1) smaller 
than the distance from solution 1. 

5.2.4. Distance with OflllJ and Na I velocities 

In Figure [9] we plot the 0[III] A5007 and Na I A5890 
RVZI values from Table [1] and again, for comparison, the 
Ha BVZI values and the radio velocity curves for the 
distance of 3.96 Mpc (solution 1). The scatter in the 
OIII RVZI values is relatively large, extending over more 
than 2000 km s'^ or ~20%. The scatter of the Na I 
RVZI values is smaller. Almost all of the velocity values 
from these two spectral lines are clearly smaller than the 
corresponding values from the Ha line, on average by 
12%. Correspondingly, the distance estimates (solutions 
7, 8) are also 12% smaller than the estimate of solution 
1. However, we think that these estimates do not reflect 
the real distance of M81 but rather the smaller maximum 
radii and velocities of the observed 0[III] and Na I line 
gas in the ejecta. 

5.3. Systematic errors 

The systematic errors of our distance estimate mostly 
depend on how reliably the radio and optical velocities 

Strictly speaking, combining the uncertainty in quadrature 
as we did gives only approximately the statistical uncertainty of 
the difference, since the data are overlapping and therefore not 
statistically independent. 



can be equated. In this context, the degree of isotropy of 
the expansion and the spatial relation between the radio 
and optical emission regions are important factors. We 
have identifled seven items concerning systematic uncer- 
tainties and elaborate on them below. 

1. Large scale anisotropy of the radio expansion: 

The degree of isotropy of the expansion of the outer 
radio shell, at least in the plane of the sky, can be deter- 
mined most directly by measuring the angular expansion 
of individual segments of the supernova from a fixed ref- 
erence point in the frame of the host galaxy. Such mea- 
surement was made for SN 1993 J relative to M81*, the 
core of the nuclear region of the galaxy, and resulted in a 
limit on anisotropic motion in any direction in the plane 
of the sky of 5.5% (Paper I). More stringent limits can be 
set on the ellipticity of the supernova's projection. For 
our composite image in Figure [l] an ellipse fitted to the 
20% contour is circular even to within 1%. We take the 
mean of these two anisotropy measurements of 3% for 
the 1(7 contribution to the systematic error. 

2. Small scale anisotropy of the radio expansion: 

As we reported in Paper HI, the composite image 
in Figure [T] shows an apparent small protrusion to the 
southwest. Such protrusions may result from Rayleigh- 
Taylor instabilities of the contact surface, which causes 
vortices through the region of the shocked CSM between 
and Oo- Such vortices are expected to show enhanced 
polarization. No significant polarization has however 
been found in several of the images we analyzed, with 
a 3(7 limit of 9% of the image peak averaged over the 
beam size. Further, the significance of the protrusion 
itself is marginal. The noise-corrected 1-a upper limit 
on the rms variation of the radius of the 20% contour 
is 3% (Paper III). Since our radio shell expansion veloc- 
ity is that of a fit spherical shell model and is therefore 
azimuthally averaged, the effect of the occurrence of pro- 
trusion would be absorbed in the statistical error of the 
expansion and does therefore not have to be considered 
as a separate element in the error budget. It needs how- 
ever to be considered again under item 4 below. 

3. Large scale anisotropy of the expansion of the optical 
line emitting and absorbing ejecta: 

In view of the almost perfect circularity of the outer 
edge of the radio shell it would be surprising if the geom- 
etry of the outer layers of the ejecta were notably differ- 
ent. By contrast, the significant linear polarization that 
was discovered in the early optical spectra and contin- 
uum indicates that some asymmetries were likely present 
in the ejecta up to a few weeks after shock breakout. 
These asymmetries were successfully modelled with dif- 
ferent geometries by Hoflich (1995), Hoflich et al. (1996), 
and Tran et al. (1997) including one with a spherical 
outer ejecta envelope and an off-center source (Hoflich 
1995), such as may be expected since the progenitor of 
SN 1993 J was a member of a binary system (Maund et 
al. 2004, see also Podsiadlowski et al. 1993). Tran et al. 
(1997) also discussed a model where the ejecta interacted 
with a clumpy and anisotropically distributed CSM. 

In any case, in the plane of the sky, the ejecta envelope 
must have been highly circular and expanding almost 
isotropically to be consistent with the near circularity of 
the radio shell from a few weeks (Bartel et al. 1994) to 
several years (Paper HI) after shock breakout. 



8 



Bart el et al. 



What is the evidence for sphericity and almost 
isotropic expansion in three dimensions? First, the good 
match between the decrease of the BVZI (abs, em) and 
the deceleration of the outer surface of the radio shell. 
Further, the almost exactly symmetrical and box-like 
shape of the Ha emission line after one year, as first 
noted by Filippenko et al. (1994). In particular, for the 
11 epochs from t — 523 to 2454 d, the mean of the ra- 
tios of the BVZI to RVZI values is -1.024 and therefore 
the BVZI and RVZI values are almost equal in magni- 
tude, which is also reflected in the small differences of 
the distance estimates of not more than 0.5 Mpc or 1.3% 
relative to 3.96 Mpc (Table[3l solutions 3, 4 relative to 1). 
In fact the late Ha line profile, a,t t > 1000 d, can be well 
fit by a spherical shell of constant emissivity (Fransson 
et al. 2005) with a shell thickness of 30% of the shell's 
outer radius. We adopt an uncertainty of our distance 
estimate due to these asymmetries of 0.05 Mpc or 1.3%. 

4. Small scale anisotropy of the expansion of the optical 
line absorbing and emitting ejecta: 

The development of Rayleigh- Taylor fingers at the con- 
tact surface certainly leads to small scale anisotropics of 
the expansion of the ejecta, but their contribution to the 
distance uncertainty is minimal. First, Rayleigh- Taylor 
fingers likely change with time. However, our distance es- 
timates at early and late times are only different from the 
estimate of 3.96 Mpc (solution 1 in Table[3]) by up to 0.05 
Mpc or 1.3%. Second, any anisotropy should be reflected 
in the differences of the BVZI and RVZI. However, the 
distance solutions with the Ha BVZI em. and RVZI val- 
ues are only different from the estimate of 3.96 Mpc by 
0.05 Mpc or again 1.3%. We therefore conclude that any 
contribution to the distance error budget due to small 
scale anisotropics is negligible or already included in the 
contribution from the large scale anisotropics. 

5. The effect of possible biases of the angular velocity fits 
In Paper II we discussed five sources of possible sys- 
tematic errors on the determinations of the angular outer 
radius of the radio shell: i) azimuthal modulation of the 
brightness along the ridge, ii) absorption in the radio 
shell center, iii) thickness of the radio shell, iv) radial 
modulation of the shell profile, and v) deviations of cir- 
cular symmetry of the shell. We estimated that the total 
resulting error of the values of 9o is less than 5% for 
the early epochs and decreasing for the later epochs and 
largely included in the errors of 6o- Here we estimate that 
the resulting error of 9o is < 1.7% for t = 30 to ~ 300 d 
and negligible for later epochs. The corresponding con- 
tribution to the error of the distance is already largely 
included in the statistical contribution, since the errors 
of the radio velocities extrapolated to the early times of 
the Ha abs. values are much larger than 1.7%. Any re- 
maining contribution is almost certainly not larger than 
half the difference between the distance solutions with 
Ha BVZI abs. values at early times and BVZI em. val- 
ues at late times. We take an uncertainty of our distance 
estimate due to these effects of 0.05 Mpc or 1.3%. 

6. The possibility of a prolate or oblate geometry for the 
supernova 

We considered the possibility that the radio shell was 
in fact prolate or oblate, but fortuitously aligned so that 
its projection in the plane of the sky was almost circu- 
lar, in other words aligned in such a way that the long 



(prolate case) or short (oblate case) axis was pointing 
toward us at a small angle. The radial expansion veloc- 
ity would then be larger (prolate case) or smaller (oblate 
case) than the transverse one, but the observed radial 
and transverse symmetry properties would be retained. 
Consider, for example, a prolate spheroid, with the long 
axis 10% larger than its other two axes. A numerical 
calculation shows that such a spheroid, randomly ori- 
ented, has only a 7% chance of being aligned so that its 
projection on the plane of the sky is circular to within 
1.4%, as is observed for SN 1993J (Paper I, see also Pa- 
per III). In this case, we would overestimate the distance 
by ~9% since the radial velocities would in fact be larger 
than the tangential ones. The error in the distance and 
the odds against fortuitous alignment are both roughly 
proportional to the deviation of the long axis from the 
others. The chances for the alignment of oblate objects 
are similarly small, but the distance error would have the 
opposite sign. Since such alignments are unlikely, we do 
not consider them further. 

7. Spatial relation between the optical and radio regions: 

The spatial relation between the radio and optical 
emission regions which is conceptually shown in Figure [2] 
is empirically fairly well supported by the correspon- 
dence between our measurements of the shell thickness 
and that predicted by analytical computations (Chevalier 
& Fransson 1994) and hydrodynamic simulations (Mio- 
duszevski et al. 2001). We think that the largest uncer- 
tainty of the spatial relation is linked to the unknown 
growth of the Rayleigh- Taylor fingers with time and the 
corresponding slight differences in the deceleration of the 
Ha line absorbing and emitting gas, the contact surface 
and the outer surface of the radio shell or the forward 
shock. While the Ha line absorbing gas (early times) is 
less decelerated within la, the Ha line emitting gas (late 
times) is less decelerated within Aa (solutions 2 and 3, 
respectively, in Table [2|) . It appears that the Ha emit- 
ting gas at its measured maximum velocity as given in 
Table [1] is continuously expanding further into the space 
between the contact surface and the forward shock front. 
This is a sign that Rayleigh- Taylor fingers are progres- 
sively eating further into the shocked CSM. 

The range of the possible extent of the Rayleigh- Taylor 
fingers determines the distance error contribution due 
to this item, provided that the gas in the Rayleigh- 
Taylor fingers remains in the range of the temperature 
for Balmer line absorption and emission in the first place. 

We determined the error contribution through bound- 
ary conditions. The upper limit was computed by assum- 
ing that at the earliest time the hydrogen gas has not yet 
expanded into the shocked CSM through Rayleigh- Taylor 
fingers but rather extends just to the contact surface. We 
solved for the distance by fitting the angular velocity of 
the contact surface to the BVZI abs. values from t = 7 to 
10 d. We get a value for the distance of 4.16 ± 0.12 Mpc. 

The lower limit was computed by prohibiting the hy- 
drogen gas velocity from exceeding the forward shock 
velocity at any time, since we have not yet seen any sig- 
nificant protrusions in the radio shell images. Because 
of the deceleration of the hydrogen gas being slightly 
weaker at late times than that of the forward shock, to 
be conservative, the computation was done for the latest 
times. We solved for the distance by fitting the angular 
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velocity of the outer radio shell surface to the Ha and 
H/3 BVZI em. values from t = 1766 to 2454 d. We get a 
value for the distance of 3.67 ± 0.08 Mpc. We therefore 
adopt a (symmetric) error of 0.25 Mpc or 6.3% as the la 
contribution to the distance error budget. 

In case of the upper- limit distance oi D = 4.16 Mpc, 
the Ha BVZI would not exceed the velocity of the contact 
surface for the first ~10 d and only marginally if at all 
exceed it for up to t ~700 d. Only from then onward 
would the Ha BVZI, and from t - 2000 d the H/3 BVZI, 
significantly exceed the (mean) contact surface velocity, 
indicative of the development of Rayleigh- Taylor fingers 
into the shocked GSM. 

In case of the lower-limit distance of 13 = 3.67 Mpc, 
the Ha absorption gas would expand with 97 ± 6% of 
the velocity of the forward shock already a few days 
after shock breakout, which would be indicative of the 
Rayleigh- Taylor fingers having reached the forward shock 
almost from the start. Also at later times, years after 
shock breakout, the Ha emission gas would keep expand- 
ing, within the errors, with the velocity of the forward 
shock. The H/3 emission gas would in fact have exceeded 
the velocity of the (mean) forward shock by as much as 
1.9(7, which is unlikely. 

5.4. The distance solution with statistical and 
systematic errors combined. 

For the statistical standard error of the distance we 
take the uncertainty from solution 1 in Table [3l For the 
systematic errors we consider our elaborations above and 
the range of solutions listed in Table [3l The largest con- 
tribution is the uncertainty of relating the locations of 
the optical and radio emission, which is largely related 
to the uncertainty of the growth of the Rayleigh- Taylor 
fingers over time. The second largest contribution comes 
from our estimates of, or limits on, anisotropies of trans- 
verse expansion. The other contributions are related to 
anisotropies of radial expansion and to changes of the dis- 
tance estimate as a function of time. They are minor in 
comparison. We list all contributions to the error budget 
in Table ID We add these latter contributions to the sta- 
tistical error in quadrature and get a combined standard 
error of 0.29 Mpc or 7.3%. Our final value of the distance 
to SN 1993J and its host galaxy M81 is therefore: 

D = 3.96 ± 0.05(stat.) ± 0.29(syst.) Mpc 

or 

£> = 3.96 ± 0.29 Mpc. 

We plot the Ha and H/5 BVZI values with the radio 
velocity curves for our upper (4.25 Mpc) and lower (3.67 
Mpc) la distance limits in Figures \W\ and [Til respec- 
tively. 

6. DISCUSSION 

The combination of our VLBI observations of SN 1993 J 
with optical spectral line observations from the time of 
explosion to several years thereafter provides a unique 
opportunity to study the spatial and dynamical relation 
between the radio shell and the optical line emitting gas 
in the context of the circumstellar interaction between 
the supernova ejecta and the surrounding hydrogen gas 
and to determine the distance to the host galaxy M81 



geometrically with the Expanding Shock Front Method 
(ESM). 

The multifrequency VLBI observations, phase- 
referenced to the core of the host galaxy, have 
significantly advanced our knowledge of the evolution 
of supernova radio shells. In our Galaxy, radio shells of 
supernovae have been observed over at most ^ 10% of 
their age. SN 1993J has been observed essentially over 
100% of its age. In the first paper of this series, and of 
relevance to this discussion, we reported the position of 
the explosion center with an accuracy of about 160 AU 
in the galactic reference frame, and determined a 5.5% 
upper limit on any anisotropic expansion on the plane of 
the sky (Paper I). In the second paper of the series we 
consistently determined the rate of the expansion of the 
supernova from the explosion center throughout most 
of the supernova's lifetime, and determined changes of 
the rate, both with high accuracy (Paper II). In the 
third paper we presented the series of images of the 
expanding radio shell and investigated the structure 
changes and the emission profile of the shell (Paper 
HI). In this fourth paper we combine our radio results 
on the speed of the expansion, the limit on anisotropy, 
and on the emission profile of the radio shell, with 
computations of the location of the contact surface, and 
optical observations of the width and shape of spectral 
lines. 

These results are important for discussions concerning 
1) optical line-emitting gas and its relation to the radio 
shell and Rayleigh- Taylor instabilities, 2) possible mis- 
interpretations of some particular Ha and He I lines re- 
ported in the literature, and 3) the distance to the host 
galaxy M81 and the Hubble constant. We will discuss 
each of the aspects in turn. 

6.1. The optical line-emitting gas and its relation to the 

radio shell 

A few aspects of the results mentioned in the previous 
section deserve further discussion. First, the RVZI and 
RVHI Ha values for i = 16 d to 500 d indicate a mod- 
erately sloped red side, with the RVZI values well below 
the velocity of the reverse shock and with the RVHI val- 
ues showing no indication of deceleration (compare Fig- 
ures [6] and [7]) . Second, for t >500 d, the Ha profile is 
very symmetric, with steeply sloped blue and red sides. 
Third, the distance estimates with the Ha abs. at early 
times and Ha em. at late times are equal to within 0.09 
Mpc or 2%. And fourth, at late times the radio shell is 
slightly more decelerated than the line emitting gas. 

The first point may indicate that at early times the red 
side of the Ha line profile is affected significantly by the 
light passing through the interior of the supernova and 
that with time and decreasing density in the interior, 
this effect on the line profile diminished. An alternative 
interpretation, that a couple of weeks after shock break- 
out, the Ha emitting gas was spread irregularly through 
a much larger fraction of the ejecta on the far side than 
on the near side of the expanding supernova, is less at- 
tractive in view of the box-like emission profile seen later 
on. The Ha line emitting gas would have to develop from 
a one-sided irregular geometry to a spherical shell with 
fairly well defined boundaries, which is unlikely. 

The second point indicates that the Ha line emitting 
region is very symmetric along the line of sight, at least 
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for t ^500 d. Indeed the high degree of symmetry, with 
the magnitudes of the maximum vefocities in the two 
radial directions being equal to within 2.4%, is similar to 
the degree of circular symmetry in the plane of the sky, 
to within 1.4% of the radio shell's outer edge, making it 
highly likely that the radio shell as well as the Ha line 
emitting gas shell are spherically symmetric. 

The third point indicates that, assuming that the max- 
imum observed velocities are indeed those of the tip of 
Rayleigh- Taylor fingers, any evolution of such fingers 
from i = 7 d to ~3000 d is, on average, rather con- 
strained. No matter whether we used a) 9o, b) 9c, or c) 
to solve for the distance, implicitly assuming that these 
fingers a) reached the forward shock, b) were largely con- 
strained to the contact surface, or c) reached to half way 
in between the contact surface and the forward shock, the 
distance estimates for t = 523 to 2454 d were larger by 
only a) 0.06±0.09, b) 0.23±0.12, and c) 0.09±0.08 Mpc 
than those for t = 7 to 19 d. 

The somewhat larger discrepancy of 0.23±0.12 Mpc 
when using 9c, could have several interpretations. It 
could be taken as an argument against the respective dis- 
tance solution of 4.33 Mpc which is anyway larger than 
our 1(7 upper limit of the distance. It could also be taken 
as information on the shell thickness at early times, the 
evolution of the Rayleigh- Taylor fingers at late times, or 
the separation of the contact surface from the forward 
shock, since using dc for the fit is, in comparison to us- 
ing 9x or 9o more dependent on these parameters. The 
shell thickness would have to be slightly larger than as- 
sumed. More precisely, for i = 7 to 19 d, rfs/rrs would 
have to be 1.22±0.03, rather than 1.16 as assumed, in 
order to give a better fit, that is the same distance for 
t = 7 to 19 d as for t = 523 to 2454 d. As an alter- 
native, the Rayleigh- Taylor fingers would have to have 
only slightly expanded into the shocked CSM toward the 
forward shock at t 1000 to 3000 d. (see Figure [TU)) . 
Lastly, it is possible that at late times the contact sur- 
face itself or the base of the Rayleigh- Taylor fingers has 
moved closer to the forward shock than assumed. Tak- 
ing into consideration any of these alternatives would 
improve the quality of the fit to that obtained using 9^ 
or 9o. In other words, there is in fact marginal evidence 
that one or more of these alternatives actually apply. 

The fourth point indicates that the comparison of the 
radio and optical deceleration may provide a sensitive 
way of probing the evolution of the Rayleigh- Taylor fin- 
gers for particular time intervals. While the radio and 
optical decelerations for the total time and at early times 
are equal within the combined standard errors, the ra- 
dio deceleration is slightly but significantly larger than 
the optical deceleration for t > 520 d. This difference 
may indeed be indicative of the Rayleigh- Taylor fingers 
eating into the shocked CSM. The 0[III] and Na I line 
emitting gas is also significantly decelerated, slightly less 
than the Ha line emitting gas, but within the errors of 
the former. Its velocity is 12% smaller and may not be as 
much influenced by the reverse shock as the Ha gas. A 
smaller deceleration, if any at all, is therefore expected. 

6.2. The distance to M81 

The nearby spiral galaxy M81 is considered an impor- 
tant one for the determination of the extragalactic dis- 



tance scale and the Hubble constant, Hq, since it is used 
as a calibrator for various methods of determining dis- 
tances and also since it can be used for a comparison 
of the various methods. How does our ESM distance 
compare with other distance determinations to M81 or 
SN 1993J? In Table[5]we compare our determination with 
those obtained with methods that depend on an abso- 
lute calibration scheme and those obtained with meth- 
ods that, like ESM, give the distance directly. These 
latter methods are the Expanding Photosphere Method, 
or EPM (Kirshner & Kwan 1974; Eastman & Kirsh- 
ner 1989; Eastman, Schmidt, & Kirshner 1996) which 
is a variation of the Baade-Wesselink method first sug- 
gested by Baade (1926) for variable stars and later mod- 
ified and applied by Wesselink (1946), and the simi- 
lar Spectral-Fitting Expanding Atmosphere Method, or 
SEAM (Baron et al. 1995). Both methods derive the 
distance by assuming spherical symmetry and, in effect, 
combining an estimate of the angular radius with the 
linear radius of the photosphere or the line-forming re- 
gion. The former further assumes that the spectrum of 
the supernova is approximately that of a blackbody and 
computes the angular radius via spectral photometry and 
the linear radius from the velocity of spectral lines and an 
estimated time of shock breakout. The latter uses non- 
LTE radiative transfer codes and fits observed spectra 
with synthetic spectra to determine the spectral energy 
distribution and the angular and linear radii of the su- 
pernova at any given time. 

The difficulty in deriving distances with EPM is that 
a velocity needs to be derived from the Doppler shift of 
those lines that are assumed to best match the photo- 
spheric velocity, and that a correction factor needs to be 
computed to account for the differences between the su- 
pernova spectral energy distribution and that of a black- 
body. Recently, Dessart & Hillier (2005) reported on new 
computations of these correction factors that would in- 
crease the EPM distances by Eastman et al. (1996) by 10 
to 20%. A direct comparison between the EPM distance 
for SN 1999em of 7.5 ± 0.5 Mpc (Hamuy et al. 2001), 
8.2 ± 0.6 Mpc (Leonard et al. 2002), and 7.83 Mpc 
(Elmhamdi et al. 2003), with the Cepheid distance for 
the host galaxy, NGC 1637, of 11.7 ± 1.0 Mpc (Leonard 
et al. 2003) shows the range of discrepancies. 

The difficulty with SEAM is to compute radiation 
transport in a rapidly expanding supernova atmosphere 
so that the resulting synthetic spectra accurately match 
the observed ones. For example, the distance to 
SN 1999em derived with SEAM was reported to be 
12.5 ± 1.8 Mpc (Baron et al. 2004), 7% larger than, 
but within the errors of, the Cepheid distance. 

In comparison to at least EPM, ESM has in principle 
more observable parameters and fewer systematic uncer- 
tainties. In ESM, the transverse expansion velocity is di- 
rectly measured and not dependent on dilution factors. 
Changes of the expansion velocity can be measured di- 
rectly. The distance determination is almost completely 
independent of the assumed date of shock breakout. The 
isotropy of the expansion can also be measured directly, 
at least in projection. What is similar to EPM and 
SEAM is the difficulty of relating the Doppler shift of 
the spectral lines to the transverse expansion velocity. 

Our small statistical uncertainty of only 1.3% reflects 
the good flt of the VLBI expansion curve to the maxi- 
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mum Doppler shift of the Ha hne. Our larger system- 
atic uncertainty of 7.2% includes a 6.3% standard error 
due to the uncertainty of combining optical and radio 
velocities, and, omitted from published EPM or SEAM 
distance determinations, a 3% uncertainty from our mea- 
surements on the anisotropy of the transverse expansion 
of the shock fronts. 

Our value of 3.96±0.29 Mpc is in the upper half of pre- 
viously determined distances with ESM and SEAM, but 
consistent with them. The Cepheid distance estimate 
of 3.93 ± 0.26 Mpc by Huterer et al. (1995) is virtually 
the same as our ESM distance estimate. In this con- 
text it is worth noting that recent distance estimates for 
the galaxy M82 and the dwarf elliptical galaxies, F8D1 
and BK5N have also been similar to our distance esti- 
mate for M81. The galaxies are all in the same group, 
and a HI map shows evidence that M81 and M82 have 
tidally interacted recently. They have a projected sep- 
aration of only ~ 40 kpc. Using the tip of the red 
giant branch method, Madore reported a distance of 
3.9 ± 0.4 Mpc for M82 and Caldwell et al. (1998) dis- 
tances of 3.98±0.15 Mpc for F8D1 and 3.80 ± 0.27 Mpc 
for BK5N. These distance estimates are in good agree- 
ment with our M81 distance estimate. 

In contrast, our distance estimate is 9±13% larger than 
the Cepheid distance estimate of 3.63 ± 0.34 Mpc by 
Freedman et al. (1994), In other words, if the distance to 
M81 were indeed 3.63 Mpc, which is only slightly below 
our 1(7 lower limit of 3.67 Mpc, then the Ha absorption 
between t = 7 and 19 d and the later Ha and H/3 emission 
at t ~ 2000 would have to occur right up to the forward 
shock (see Figure [TO)) . The most likely interpretation for 
this scenario would be that the Rayleigh- Taylor fingers 
developed strongly right after shock breakout, reached 
the forward shock at t ~ 10 d and and continued to 
reach the forward shock up to at least t 2500 d. We 
consider this scenario to be less likely and instead argue 
in favor of a larger distance. 

Direct and Cepheid distance determinations can also 
be compared for the nearby galaxies M33 and NGC 4258. 
Using a detached eclipsing binary, the distance to M33 
was found to be 0.96±0.05 Mpc (Bonanos et al. 2006), 
14% or 1.5 X the combined uncertainty larger than the 
Cepheid distance of 0.84±0.06 Mpc (Freedman et al. 
2001). In contrast, the distance derived with maser VLBI 
is 0.73to \l Mpc (Brunthaler et al. 2005), 13% or 0.6x 
the combined uncertainty smaller than the Cepheid dis- 
tance. For NGC 4258, using maser VLBI, the distance 
is 7.2 ± 0.5 Mpc (Herrnstein et al. 1999). This value is 
10% or Ix the combined uncertainty smaller than the 
Cepheid distance of 7.98 Mpc (systematic uncertainty 
not given, total error assumed to be 0.6 Mpc, i.e., the 
same as that given by Newman et al. 2000; compare with 
8.1±0.4 Mpc by Maoz et al. 1999 and 7.8 ± 0.6 Mpc by 
Newman et al. 2000). It appears that, at least for our 
examples, distances derived with ESM, SEAM, and by 
using a detached eclipsing binary, are larger than the 
Cepheid distances by 7 to 14% while the distances de- 
rived with maser VLBI tend to be somewhat smaller. 
Perhaps the systematic uncertainties are still underesti- 
mated. 

6.3. The Hubble Constant 



Is it possible to derive Hq from our distance estimate 
alone? Recently it has been reported that the Hub- 
ble flow can be extrapolated back to about 2 Mpc if it 
is corrected for peculiar motions caused by the large- 
scale distribution of matter in the nearby universe (e.g., 
Karachentsev et al. 2003; Sandage et al. 2006). Sev- 
eral authors give corrected flow velocities for M81 or the 
M81 group: Tonry et al. (2000) give 246 km s'^ (given 
in Freedman et al. 2001 for M81) and Sandage et al. 
(2006) give 234 km s~^. With our distance estimate, 
these flow velocities would lead to Hq of 62 ± 5 and 59 ± 4 
km s~^ Mpc~^, respectively. Karachentsev et al. (2002) 
give a low velocity, of 106 km s~^, but a rather large 
velocity of 360 km s~^ for the spatially nearby galaxy 
M82. Freedman et al. (2001) give a similarly low veloc- 
ity for M81, of 80 km s^^. Such low velocities would lead 
to unrealistically small values of Hq. Clearly, the range 
of predicted flow velocities for M81 is still too large to 
be used for a reliable estimate of Hq and our distance 
estimate for M81 can at best be taken as an indication 
against a relatively large value for Hq. 

Is it possible to derive Hq by anchoring the Cepheid 
distance scale to our distance estimate for M81? In this 
respect the Cepheid distance modulus of M81 relative to 
that of the Large Magellanic Cloud (LMC) may be of 
interest since the LMC is in many studies crucial for the 
entire Cepheid distance scale. However, different authors 
derive different relative distance moduli. For instance, 
Freedman et al.'s (1994) distance modulus for M81 of 
27.80 ± 0.20 mag, corresponding to a distance, Dmsi — 
3.63 ± 0.34 Mpc, is based on the assumption of an LMC 
modulus of 18.50 mag (Dlmc = 50.1 kpc). 

In contrast, Huterer et al. (1995) obtained, with a dif- 
ferent analysis of the same HST data, an M81 modulus 
of 27.97±0.14mag (Dmsi = 3.93±0.26 Mpc). This esti- 
mate, as mentioned earlier, is essentially identical to our 
estimate for the distance of M81. In the same analysis 
they also obtained an LMC modulus of 18.45 ± 0.10 mag 
(^'lmc = 49.0 ±2.3 kpc). Apparently, different analyses 
led to distance ratios different by 10%. 

Huterer et al.'s (1995) LMC distance modulus can 
be compared with that of 18.41 ± 0.16 mag obtained 
by Macri et al. (2006) through a determination of the 
Cepheid distance modulus of the maser-host galaxy NGC 
4258 relative to that of the LMC and who give an esti- 
mate of i?o = 74 ± 3(stat.) ± 6(syst.) km s"^ Mpc"^ 
It can also be compared with 18.39 ± 0.05 mag ob- 
tained through revised Hipparcos parallaxes of Cepheids 
by van Leeuwen et al. (2007) who also argue in favor 
of a revision of Sandage et al.'s (2006) and Freedman 
et al.'s (2001) values for the Hubble constant to Hq = 
70 ± 5 km s-i Mpc-i and i/o = 76 ± 8 km s'^ Mpc"\ 
respectively. 

Which value for Hq can we derive with our dis- 
tance estimate for M81 in view of this discussion? 
Freedman et al. (2001) obtained, on the basis of their 
Cepheid analysis, Hq = 72 ± 8 km s'^ Mpc~^ If 
the difference between their M81 distance determina- 
tion and ours is due to a systematic effect, then, on 
the basis of their Cepheid analysis, scaling would lead 
to Hq — 66 ± 11 km s~^ Mpc~^, where the uncer- 
tainty is one standard error, derived by adding errors 
in quadrature. This value can be compared with Hq — 
62.3±1.3(stat.)±5.0(syst.) km s^^ Mpc^^ (Sandage & 
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Tammann 2006). Hutcror ct al.'s (1995) Ccpheid dis- 
tance determinations for M81 and the LMC, however, 
are consistent with a ~ 10% higher value for Hq. 

In conclusion, our distance estimate for M81, which is 
somewhat larger than Freedman et al.'s (1994) estimate, 
may argue for a correspondingly smaller value of Ho. 
However, deriving a reliable value of Hq from our M81 
distance determination directly would require a more re- 
liable estimate of MBl's flow velocity than is presently 
available. Deriving such value of Hq indirectly would re- 
quire anchoring other distance scales more reliably to our 
M81 distance determination. 

7. CONCLUSIONS 
Here we list a summary of our main conclusions. 

1. SN 19993 J is the best suited supernova yet for deter- 
mining an accurate geometric distance with the expand- 
ing shock front method. ESM, by comparing velocities 
from the Doppler shift of optical lines with the VLBl 
determined velocities of the radio shell. 

2. The brightness profile of the radio shell being consis- 
tent with a sharp intensity decline at the outer rim ar- 
gues in favor of the outer surface of the radio shell being 
indeed equivalent to the expanding forward shock front. 
The thickness of the radio shell being consistent with pre- 
dictions from hydrodynamic simulations argues in favor 
of the inner surface of the radio shell being equivalent to 
the reverse shock front. 

3. The angular expansion velocity of the outer surface 
of the radio shell could be determined for the time from 
t ~ 30 to > 3000 d with an uncertainty at early times of 
<5% and at late times of <2%. 

4. The average of the thickness of the radio shell, which 
is 29% of the radius of the outer radio shell surface for 
the time t 1000 to 3000 d, coincides with the average 
of the thickness predicted by hydrodynamic simulations 
for that time interval. 

5. Changes of the velocity of the radio shell's outer sur- 
face at i ~ 350 and 2000 d, which are reflected by sharp 
minima in the X-ray lightcurve at these same times, have 
only partial correspondence in the velocities from the op- 
tical lines. At t ~ 350 d, only Ha RVZI and RVHI were 
recorded and do not show any correspondence, providing 
evidence that the observed Ha velocities at the far side of 
the expanding supernova were biased by absorption and 
do not reflect the velocity of the contact surface. At the 
less dominant change at t ~ 2000 d, only the velocities 
from the H/3 and 0[111] lines appear to be influenced. 

6. The Ha gas has with m = 0.868 ± 0.004 virtually 
the same (mean) deceleration as the outer surface of the 
radio shell (forward shock, m = 0.870 ± 0.005) and the 
surface midpoint between the outer surface and the con- 
tact siuface (m = 0.864 ± 0.008) for the total time be- 
tween t = 7 and ~2500 d. The decelerations are also 
equal within the errors for the early time interval be- 
tween f. = 7 and ^300 d. For the late time interval from 
t - 500 to -2500 d the Ha gas is with m = 0.890 ±0.022 
slightly less decelerated than any of the surfaces of the 
radio shell. 

7. It is possible that the difference in the deceleration is 
due to the development of Ray leigh- Taylor fingers. 



8. The 0[I11] and Na 1 gas is also decelerated, although 
only with marginal significance. 

9. Using ESM by combining the BVZI values of the 
Ha line with the angular velocity half way between that 
of the contact surface and that of the radio shell's outer 
surface gives a geometric distance estimate of D = 3.96± 
0.29 Mpc. The standard error combines in quadrature 
a statistical contribution of 0.05 Mpc and a systematic 
contribution of 0.29 Mpc. 

10. The distance estimate depends by not more than 
0.06 Mpc or 1.5% on whether the estimate is obtained 
with only a) early Ha BVZI data from < = 7 to 19 d, 
b) late Ha BVZI data from t - 500 to ~ 2500 d, c) Ha 
RVZI data, or d) H/3 BVZI data. 

11. The largest contribution to the error comes from the 
uncertainty of relating the optical velocities to the radio 

velocities. 

12. The l(j upper limit of 4.25 Mpc corresponds to the 
Ha absorption occurring at the earliest time (t — 10 d) 
only up to the reverse shock. The la lower limit of 3.67 
Mpc corresponds to the Ha absorption or emission and 
H/3 emission occurring as far out as the forward shock 
which would mean that Rayleigh- Taylor fingers have de- 
veloped and stretched out to the forward shock within 
the errors, from essentially t ^ 10 to — 2000 d. 

13. The distance estimate is 9 ± 13%, or 0.7 x the 
combined errors, higher than the Cepheid distance of 
3.63±0.34 Mpc by Freedman et al. (1994). 

14. If the distance to M81 were indeed 3.63 Mpc, then 
the Ha absorption at t ^ 10 d and the Ha and H/? emis- 
sion at t ~ 2000 d would have to originate as far out 
as the forward shock. The most reasonable interpreta- 
tion would be that Rayleigh- Taylor fingers had developed 
right after shock breakout and stretched out as far as the 
forward shock. We think that this scenario is less likely 
and that the distance has to be larger. 

15. Our direct distance estimate for M81 may argue for 
a value of Hq somewhat smaller than that of Freedman 
et al. (2001). However, to derive a more reliable value 
of Hq from our M81 distance determination, a more re- 
liable Hubble flow velocity estimate or a reliable way of 
anchoring other distance scales to the distance of M81 
would be needed. 
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APPENDIX 

POSSIBLE MISINTERPRETATIONS OF SOME PARTICULAR Ha AND HE I LINES BY OTHERS 

Matheson et al. (2000b) reported a couple of Ha BVZI values of 16,600 and 16,100 km at t = 433 and 473 d, 
respectively, and He I A5876 BVZI values of ~16,000 km s~^between t = 553 and 2454 d that appear puzzling in 
the context of the other optical velocities and the radio shell velocities. Each of the values is clearly larger than the 
velocity of the outer surface of the radio shell, at late times even twice as large. The couple of Ha BVZI values are 
almost certainly overestimates due to blending of the Ha line with the 0[I] line. On first sight it is intriguing that 
the apparent large decrease in the Ha velocity from 16,600 and 16,100 km s^^ a.t t ~ 433 and 473 d to 10,400 and 
10,000 km s~^ at t = 523 and 533 d, respectively, approximately coincides with the large deceleration of the radio shell 
which in turn coincides with an increase in the X-ray flux (Zimmermann & Aschenbach, 2003). However, we think 
that the diminishing effect of the blending of the Ha line at the time of the large radio deceleration is coincidental. 
The He I values were derived from a decomposition of the line profile. This decomposition may be questionable in 
light of the values being so large. We did not list these Ha and He I velocity values in Table [T] 
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TABLE 1 

Expansion velocities of SN 1993J from optical lines 



Age-* 


Ha A6563 


Ha A6563 


Ha A6563 


Ha A6563 


H/3 A4861 


[OHI] A5007 


Na I A5890 


Reference*^ 




Jj V ZjL 


Jj V iTi 


ixV ZjI 


rlV xll 


Jr> V ZjI 


IXV ZjI 


IXV ZjI 




(d) 


(km s-1) 


(km s-l) 


(km s-i) 


(km s"^) 


(km s-i) 


(km s-i) 


(km s-i) 




7 


-19,200±500'* 














1 


8 


-19,000±500'* 














1 


9 


-18,900±500'* 














1 


10 


-18,000±500'' 














1, 2, 4 


16 


-18,000±500'' 




14,000ib500 


6,000ib500 








2,4 


17 


-18,000±500'' 




14,000ib500 


6,000ib500 








2, 4 


18 


-18,000±500f 




14,000ib500 


6,000ib500 








2, 4 


19 


— 18,000±500'^ 




14,000di500 


6,000di500 








3 4 


123 






12,600ib500 


5,500ib500 








4 


171 






ll,420di400S 








5 


182 






11 min-i-^nn'^ 










4 


236 
















5 


zoo 



































4 


298 






11 c;nn-i-'^nn 


7 zinn-i-'^nn 








4 


299 
















5 


315 






11 c;nn-i-'^nn 










4 


367 






1 1 zLfin-i-4nns 










5 


355 








7 (sno-i-'^nn 








4 


387 








8 200+^00 

(J J iiUW J_!JUVJ 








4 


433 








a 200+50(1 








4 


473 






10,700±500 


8,800±500 








4 


523 


-10,900±500 


-9, 100±500 


10,400±500 


8,400±500 








4 


553 


-10,700±500 


9, 100±500 


10,500±500 


8,200±500 


-10,500±1,050 


9,100±910 


8,400±840 


4 


670 


-10,500±500 


-8, 800±500 


10,500±500 


8,200±500 


-9, 700± 970 


9,600±960 


9,200±920 


4 


881 


-10,400±500 


-8, 600±500 


10,000±500 


8,200±500 


-9, 500± 950 


10,000±1,000 


8,700±870 


4 


976 


-10, 200±500 


-8, 700±500 


9,600±500 


7,900±500 


-9, 500± 950 


8,000±800 


8,900±890 


4 


1766 


-9, 600±500 


-8, 300±500 


9,300±500 


7,700±500 


-9, 700± 970 


7,600±1,520 


8,500±850 


4 


2028 


-9, 000±500 


-7, 600±500 


9,500±500 


7,300±500 


-10,400±1,040 


9,800±980 


8,400±840 


4 


2069 


-9, 200±500 


-7, 600±500 


9,100±500 


7,300±500 


-12,000±1,200 


7,500±870 


8,700±870 


4 


2115 


-9, 300±500 


-7, 600±500 


9,000±500 


7,300±500 


(-14,900±2,980) 


7,200±1,440 


8,300±830 


4 


2176 


-9, 300±500 


-7, 800±500 


8,700±500 


7,400±500 


(-11,800±2,360) 


6,700±1,340 


7,700±770 


4 


2454 


-9, 300±500 


-8, 200±500 


9,200±500 


7,200±500 


(-13,400±2,680) 


8,400± 1,680 


7,500±750 


4 



^ Time since assumed explosion date of 1993 Mar. 27.5 UT. Note, that this time is 0.5 d earUer than the time we assumed in general throughout 
this paper and Papers I, II, and III. The difference is not significant in any of our computations.*^ (1) Finn et al. 1995; (2) Lewis et al. 1994; (3) 
Trammell et al. 1993; (4) Matheson et al. 2000a; (5) Patat et al. 1995.*^ BVZI: maximum velocity on the blue edge of the line; BVHI: velocity at 
half of the maximum intensity on the blue side of the line: RVZI: maximum velocity at zero intensity on the red edge of the line; RVHI: velocity on 
the red side of the line at half of the maximum intensity. The lines are emission lines except on days 7 to 19, see (d, c). We determined the velocity 
values from the published line profiles, except where indicated otherwise, see (f). The uncertainties are estimated standard errors, given by the 
uncertainty of determining the baseline of the profiles, the values at the baseline, the maximum of the profile and the values at the half maximum. 
For more information, see text.*^ The velocity values and their errors were taken as listed in the literature. The values in parentheses have extra 
large errors and were not considered further in our analysis.'^ Maximum velocities obtained by adding 1400 km s""'" to the reported velocities for 
the minimum of the absorption trough of 17,800, 17,600, and 16,800 km s""'" (Finn et al. 1995). The offset of 1400 km s""'" was obtained from the 
difference of the reported velocity for the minima of the troughs of 16,600 km s""'" (Finn et al. 1995) and the measured maximum velocity of 18,000 
km s~^ from a published line profile (Lewis et al. 1994).^ Maximum expansion velocity from the blue edge of the absorption trough of the Ha 
line.^ Values listed by Patat et al. 1995. Profiles were not plotted for all these epochs.^ Only the RVZI value is given. The RVHI value is omitted 
because blending of the 0[I] line made identifying the maximum of the Ha profile uncertain by a large margin. 
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TABLE 2 

The deceleration parameter for the radio and optical data 



4t 
TT 


Radio 


range 




Optical 






Lines'' 








t_/ LtH-jOjI v_l.t,v_.l:,l • 




t(d) 


t(d) 


v(kja 






mo 


TTlc 


nix 


"lopt 


1 


30 - 


2,432 


7 - 


2,454 


19,200 - 


- 9,000 


Ha abs., em. (BVZI) 


0.870 ± 0.005 


0.855 ± 0.008 


0.864 ± 0.008 


0.867 ± 0.004 


2 


30 - 


- 306 


7 - 


- 19 


19,200 - 


18,000 


Ha abs. (BVZI) 


0.919 ±0.019 


0.904 ± 0.020 


0.912 ± 0.020 


0.936 ±0.023 


3 


520 - 


2,432 


523 


- 2,454 


10,900 - 


- 9,000 


Ha em. (BVZI) 


0.798 ±0.007 


0.782 ± 0.009 


0.791 ± 0.009 


0.890 ± 0.022 


4 


520 - 


2,432 


523 


- 2,454 


10,700 - 


- 8,700 


Ha em. (RVZI) 


0.798 ± 0.007 


0.782 ±0.009 


0.791 ± 0.009 


0.897 ±0.023 


5 


520 - 


2,432 


553 


- 2,454 


12,000 - 


- 9,500 


H/3 em. (BVZI) 


0.798 ±0.007 


0.782 ±0.009 


0.791 ± 0.009 


1.062 ±0.063 


6 


30 - 


2,432 


7 - 


2,454 


19,200 - 


- 9,000 


all of # 1 - 5 


0.870 ±0.005 


0.855 ± 0.008 


0.864 ± 0.008 


0.865 ± 0.003 


7 


520 - 


2,432 


553 


- 2,069 


10,000 


- 7,500 


OpH] em. (RVZI) 


0.798 ± 0.007 


0.782 ±0.009 


0.791 ± 0.009 


0.908 ± 0.066 


8 


520 - 


2,432 


553 


- 2,454 


9,200 - 


7,500 


Na I em. (RVZI) 


0.798 ± 0.007 


0.782 ±0.009 


0.791 ± 0.009 


0.929 ±0.048 



^ As in Tablc[3^ The deceleration parameter, m, with r oc and v oc mi"^^^, for the outer radius of the radio shell (mo), the radius of the 
contact surface (m-c), the radius midpoint between them (mx), and the velocity of the ejecta gas with optical line absorption and emission (rriopt)- 
Sec ^ |4.3| for more information. The solutions of rrio have a Chi-square per degree of freedom, 5 close to unity. The errors are adjusted for — 1- 
The errors in mc and mx are taken from those of mo added in quadrature to 0.005 assumed to be the error from the theoretical models. The 
solutions of mopt have x^ between 0.1 and 0.7 for the hydrogen gas, 1.5 for the 0[III] gas and 0.3 for the Na I gas. The errors are not scaled for 
X^ — 1, since we think that our measurements of the velocities from the line profiles are somewhat correlated. For solution 1 and 6, we used data 
on t — 30 d and between t — 520 and 2,432 d to best match the time range in which optical data were used. The solution was forced to fit the data 
point at i — 30 d exactly to give appropriate weight to this early point. 



TABLE 3 

Distance estimates from combined radio and optical velocity data 



Solution 


Optical data range^ 


Lines'' 


Fit distance'^ 


No. of data points 


2d 


# 


t 


V 




D 








(d) 


(km s~^) 




(Mpc) 






1 


7 - 2,454 


19,200 - 8,800 


Ha abs., em. (BVZI) 


3.96 ± 0.05 


19 


0.6 


2 


7 - 19 


19,200 - 18,000 


Ha abs. (BVZI) 


3.92 ± 0.06 


8 


0.2 


3 


523 - 2,454 


11,100 - 8,800 


Ha em. (BVZI) 


4.01 ± 0.07 


11 


0.9 


4 


523 - 2,454 


10,800 - 8,600 


Ha em. (RVZI) 


3.91 ± 0.07 


11 


0.9 


5 


553 - 2,454 


12,000 - 9,500 


H/3 em. (BVZI) 


4.02 ± 0.17 


7 


1.6 


6 


10 - 2,454 


18,000 - 8,600 


all of # 1 - 5 


3.95 ± 0.04<= 


37 


0.9 


7 


553 - 2,069 


10,000 - 7,500 


0[III] em. (RVZI) 


3.48 ± 0.14 


6 


1.2 


8 


553 - 2,454 


9,200 - 7,500 


Na I em. (RVZI) 


3.48 ± 0.11 


10 


0.5 



^ Time and velocity ranges of the optical data used for the fit.^ The optical lines from Table [^which were taken for the fit for specific distance 
solutions. The type of lines are Ha in absorption (abs.) and emission (em.) and H/3 in emission (em.). The velocities for absorption are BVZIs 
measured at the blue edge of the absorption trough (Lewis et al. 1994) or the blue edge of the polarization profile (Trammell et al. 1993). The 
velocities from the emission profiles are BVZIs measured at the blue edge of the profile (Matheson et al. 2000a) and RVZIs measured at the red 
edge of the profile (Patat et al. 1995).^ The distance solutions from the weighted least-squares fits of the radio angular velocities to the optical 
velocities. For the fits the time ranges of the radio and optical data were chosen to be essentially equal. For the time range 7 < i < 19 d the radio 
angular velocities were extrapolated from later values (see text). The uncertainties are statistical standard errors. Note: the errors are not scaled 
to x^ (see last column) of unity. Chi-squared per degree of freedom.^ Solution from all data used in the previous solutions combined. 



TABLE 4 
Distance error budget 



# 


Error component 


a 


Error source 




(Mpc) 


{%) 




1 


0.05 


1.3 


statistical'' 


2 


0.12 


3.0 


anisotropic transverse expansion'^ 


3 


0.05 


1.3 


anisotropic radial expansion"^ 


4 


0.05 


1.3 


difference from BVZI abs. and em."^ 


5 


0.25 


6.3 


radio-optical relation from boundary conditions' 


2-5 


0.29 


7.2 


rss^ of systematic error contributions'^ 


1-5 


0.29 


7.3 


rss of statistical and systematic error contributions 



^ All error components are standard errors.^ From solution 1 in Table [3] largely includes error due to small-scale anisotropic transverse expansion 
(see item 2. in ^ |5.3| l ."^ Largely includes error due to large-scale anisotropic transverse expansion (see item 1. in ^ |5.3l l ."^ Difference between distance 
solutions with BVZI (solution 3 in Table [3} and RVZI values (solution 4 in Table [sj. This contribution largely includes error due to large-scale 
anisotropic radial expansion (see item 3. in §[5?3j.'^ The larger of the differences of the distance solutions with BVZI abs. and BVZI em. values to 
solution 1 (solution 3 minus solution 1 in Table[3}. This contribution includes systematic errors of the angular velocity fits (see item 5. in ^ I5.3| l . 
^ The uncertainty of relating optical and radio velocities (see item 7 in ^ |5.3| l.^ Root-sum-square. 
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TABLE 5 

Recent determinations of the distance to M81 



# 


Method 


Distance (Mpc) 


Reference 


1 


ESM 


3.96±0.29 


this work 


2 


HST Cepheids 


3.63±0.34 


Freedman et al. (1994) 
Huterer et al. (1995)'' 


3 


HST Cepheids 


3.93±0.26 


4 


PN luminosities'' 


3.5 ±0.4 


Jacoby et al. (1989) 


5 


IR Tully-Fisher<= 


3.7 ±0.5 


Aaronson, Mould, & Huchra (1980) 


6 


El'.Vld 


2.6 ±0.4 


Schmidt et al. (1993)<= 


7 


EPM 


4.2 ±0.6 


Wheeler et al. (1993)'' 


8 


EPM 


3.65±1.45S 


Prabhu et al. (1995) 


9 


EPM 


3.5 ±0.2 


Clocchiatti et al. (1995)'' 


10 


SEAM' 


4.0 ±0.5 


Baron et al. (1995) 



^ Based on the Cepheid observations of Freedman et al. (1994). Planetary nebula luminosities.*^ Infrared Tully-Fisher method. Expanding 
photosphere method.^ Eastman et al. (1996) noted that a peculiar Type II supernova like SN 1993J may not be a good candidate for either 
applying EPM or assessing its accuracy.*^ A note is given in the paper that the result depends on the detailed structure of the atmosphere.^ The 
published value is 2.2-5.1 Mpc.^ The complete equation for the distance is D = ^(3.48 — 0.09<5to + 0.09(5i?o =t 0.20) Mpc where the scattering 
parameter C ^ 1) and where 5to is the difference between the actual time of shock breakout and JD=2,449,074.6 in tenths of a day and SRq is the 
difference between the actual radius of the progenitor star and 3.86x10^^ cm, in units of 10^^ cm.^ Spectral-fitting expanding atmosphere method. 




Fig. 1. — Left panel: An optical image of the spiral galaxy M81 (Sandage 1961) with a VLA contour image overlaid and with SN 1993J 
prominent as the brightest radio source after the nucleus of the galaxy. Right panel: A composite VLBI image of SN 1993J at 8.4 GHz 
from three epochs at t = 2080, 2525 and 2787 d (1998 December to 2000 November) all corrected for position shifts, scaled in flux density 
and radius to the values of the 2000 November data, and then combined (adopted from Paper III). The beam (FWHM, 0.70 mas) is shown 
at lower left. The contours are drawn at —16, 16, 32, 45.3, 64 and 90% of the peak brightness, and the background rms was 5.3% of the 
peaJt brightness. North is up and east to the left. 
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Photosphere 




Neutral ejecta, T~ 10 K 
Ionized ejecta, T~10^ K 
Reverse shock 

7 

Shocked ejecta, T~10 K 
Contact surface with Rayleigh-Taylor fingers 
Shocked CSM, T~3xlO^K 
Forward shock 
CSM, T~3xlO^K 

Radio shell model emission 




Radial brightness profile 



Model radial brightness profile 



Model radial emissivity profile 



4 

e. 



2 

9(mas) 



Fig. 2. — Top panel: The left part of the composite VLBI image of SN 1993J from Figuref^ The outer overlayed eircle shows the fit outer radius, 
^o, of the shell model with absorption in the center and indicates the expected location of the forward shock front, which is expanding into the 
circumstellar medium (CSM). The inner circle shows the fit inner radius, ^i, of that shell model and indicates the expected location of the reverse 
shock front. The forward and reverse shock fronts travel in co-moving opposite directions (see arrows) from the contact discontinuity, or contact 
surface, where the ejecta hit the CSM. The location of the contact surface is shown by the dashed circle. Fingers expanding into the shocked CSM 
due to the contact surface being Rayleigh-Taylor unstable are also shown. The smaller dashcd-dotted circle indicates the boundary between the 
colder, neutral ejecta and the hotter, ionized ejecta. Broad line emission is expected from the ionized ejecta up to the contact surface. The center of 
the circles is at the fit center position (see Paper I). Typical temperatures are also indicated. Second panel from top: The brightness profile of the 
composite image above, averaged over all p.a.'s, plotted as a function of angular radius, 9. The resolution is 0.70 mas, and the plotted uncertainties 
are the standard errors of the bin values, derived from the number of beam areas within each bin and the larger of the standard deviation within 
that bin and the rms of the noise of the background brightness. The plotted values are correlated, especially at small radii, because they are less 
than 1 beamwidth apart. Also indicated is the corresponding radial profile of the projected best fit spherical shell model with uniform emissivity and 
an absorption disk in the center, convolved to the resolution of 0.70 mas. Third panel from top: The corresponding radial profile of the projected 
unconvolved shell model with absorption. Bottom panel: The radial profile of the volume emissivity of the shell (without absorption). The shell is 
limited by 6^ and 9i. 
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Fig. 3. — The ratio 9o/di, which is a measure of the shell thickness, as a function of time. The data points are our measurements at 
times when the shell was sufficiently large for a thickness determination (from Paper II). The thin horizontal line from t = 996 d to 2996 d 
gives the mean of do/di = 1.29 for this time interval. The dotted line gives the prediction from hydrodynamic simulations. 



CO 
CO 



• 1— I 

cd 

CD 
O 

3 

< 



0.1 r 



T 1 1 1 — I — T 



1 1 — I — I — I — r 



o 


2 2 . 2 


GHz 


A 


14.8 


GHz 


• 


8.4 


GHz 


□ 


5.0 


GHz 


o 


2.3 


GHz 


o 


1.7 


GHz 




J I I I I I L 



100 1000 
Time since shock breakout (d) 



Fig. 4. — The angular outer radius, Oo, of the SN 1993J shell model as a function of time since shock brealcout. TaJsen from Paper II. 
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T — I — I — I — r-| 1 1 1 1 — I — I — I — r"| 1 1 1 ] — I — I — I — 




Time since shock breakout (d) 



Fig. 5. — The time derivatives of the angular radii of the outer {9o, upper solid curve) and inner lower solid curve) surfaces of the 
radio shell, the contact surface (9c, curve with short dashes), and the mean between the outer surface and the contact surface (9x, curve 
with dots). The thin curves with long dashes on both sides of the upper solid curve indicate the (asymmetric) standard error of 1 to 5% 
of as a function of time (see text). 
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Fig. 6. — The observed expansion velocities from the width of the Ha line from Table [T] as a function of time. The symbols represent 
the BVZI (blue velocity at zero intensity) values from the absorption trough of the line (filled square), and the BVZI (filled circle), BVHI 
(blue velocity at half intensity, filled triangle), RVZI (red velocity at zero intensity, open circle) and RVHI (red velocity at half intensity, 
open triangle) values, all from the emission line. 




Fig. 7. — As in Figure |6] but now showing only the BVZI Ha values to which DB^ (Mean) was fit to determine the distance, D. In 
addition the velocity curves, D6o, D8i, D6c, and D9^ are plotted (see Figure[31l for the best fit distance of D = 3.96 Mpc. 
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Fig. 8. — As in Figuro[7]but now also with the BVZI values of the H/3 emission line from Table [Tl 
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Fig. 9.— 



As in Figure [T] but now also with the RVZI values of the 0[III] and Na I emission lines from Table [T] 
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Fig. 10. — As in Figure[7] but now the velocity curves are plotted for a distance of 3.67 Mpc, which is the Icr lower limit of our distance 
estimate of 3.96 ± 0.29 Mpc. 




Fig. 11. — As in Figure[7] but now the velocity curves are plotted for a distance of 4.25 Mpc, which is the la upper limit of our distance 
estimate of 3.96 ± 0.29 Mpc. 



